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^.^' In this paper we study the transient surface cavity which is created by the controlled 

'"^ . impact of a disk of radius ho on a water surface at Froude numbers below 200. The 

dynamics of the transient free surface is recorded by high speed imaging and compared 
to boundary integral simulations. An excellent agreement is found between both. The 
flow surrounding the cavity is measured with high speed particle image velocimetry and 
is found to also agree perfectly with the flow field obtained from the simulations. 
We present a simple model for the radial dynamics of the cavity based on the collapse of an 
^.,^. infinite cylinder. This model accounts for the observed asymmetry of the radial dynamics 

(~| ' between the expansion and contraction phase of the cavity. It reproduces the scaling of 

Q^, the closure depth and total depth of the cavity which are both found to scale roughly 

as ex Fr ' with a weakly Froude number dependent prefactor. In addition, the model 
accurately captures the dynamics of the minimal radius of the cavity, the scaling of the 
^ ' volume Vbubbie of air entrained by the process, namely Vbubbie/hQ ex (1 + 0.26Fr^'^)Fr^ ", 

QO , and gives insight into the axial asymmetry of the pinch-off process. 

p 

f^ ' 1. Introduction 

~j^ • A spectacular example of free surface flow is the impact of an object on a liquid: 

The impact creates a splash and a transient cavity. This surface cavity then violently 
collapses under the influence of the hydrostatic pressure. At the singularity where 

►v> , the walls of the cavity collide, two powerful jets develop, one downwards and the 

}_( ' other one upwards up to several meters high, making this fast event an impressive 

5^ I scene. Research into the physics of these transient surface cavities started at the 

beginning of the twentieth century when A.M . Worthington published his famous 
work " A study of splashes" ( WorthingtonI ( 19p8[)). His photograph s revea led a wealth 



of phenomena of unanticipated complexity (jWorthington &: Cold (jlSQ?! )'). Although 



much has been contributed to the understanding of these phenomena, many o f the 



intrigu i ng questions posed by Wor thington's photographs resonate still today (jRein 
(|l993l ): lFedorchenko &: Wa^ (|2004l) '). 



All investigations since Worthington's studies entailed experiments with a freely 
falling object impacting on the free surface. To gain further insight into such impact 
events, we built a setup in which wc attach the impacting object to a linear motor. In 
this way we gain full control over the impact velocity, which now turns from a response 
observable into the key control parameter of the system. 
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The dynamics of a surface cavity are of enormous practical importance in many 
natural and indust r ial p ro cesses: Rai n drops falling onto the o c ean entrain air 
(|Oguz fc Prosperettil (|l990f) : lOguz et al\ (|l995l ): IProsperetti k Qgud (|l997t )) and it is 
this mechanism which is one of the major sinks of carbon dioxide from the atmosphere. 
Droplet impact a nd the subsequen t void collapse are also a significant source of 
underwater sound (jProsperetti et al\ (|l989l )) and a thorough understanding is therefore 
crucial in sonar research. High speed water impacts and underwater cavity formation 



Lee et 



are m or eover of relevance to military operations (Gilbarg fc AndersonI ()194i 

(|l997l ): iDuclaux et al\ (|2007h ). In the context of industrial applica tions, drop impac t 
and the subsequent void formation are crucial in pyrometallurgy i Liow et al\ ( 1996 h 
Morton et al. ( 20001) ). in t he foo d industry, and i n the c o ntext of ink-jet printin g (JLel 
(|l998l ): IChen fc BasaraT] (J20o3); Ide Jong et oil (|2006ah : Ide Jong et all pOOebh ). A 
similar series of events as in water can even be obse r yed when a s t eel ba ll impacts on 
very f i ne and soft sand (IThoroddsen fc ShenI (12001); iLohse et al\ (|2004f ): iRover et 
(|2005l ): ICaballero et al\ (|2007l) ). 



Although in some of the literature the deceleration of the impacting body was 
minimized by choosing the properties of the body such that the velocity of the im- 
actor remained roughly constant durin g the time the cavity dynamics were observed 
Glasheen fc McMahon (1996 ); Gaudetj ( 19981 )'). the velocity of the body nevertheless 
remained a response parameter set by the system. Our use of a linear motor to accurately 
control the position, velocity, and acceler ation of the im pacting o bject co nstitutes the 
key difference between our work (see also iBergmann et al. (200(3i) ; iGekle et al. (20081 )) 
and all previous literature. 



In this article, we will use observations from experiments and boundary integral simu- 
lations to construct a model which accurately describes the radial dynamics of the cavity. 
In Section [2] we present results from our controlled experiment and compare them to the 
boundary integral simulations. More specifically, in subsection 12 .31 we discuss the dynam- 
ics of the free surface and continue in subsection 12.41 with the topology and magnitude 
of the flow surrounding the cavity obtained by particle image velocimetry. 
In Section [3] we will derive a model which captures the radial dynamics of the cavity. We 
will use the model to investigate the following key characteristics of the transient surface 
cavity: First, the depth at which the pinch-off will occur is discussed in subsection 14.11 
Then, in subsection 14.21 the amount of air entrained by the cavity co llapse is studied 
The a rticle is concluded in Section [51 The results from our earlier paper iBergmann et 



( 20061 ) that are relevant to the present study are reviewed in Appendix A, together with 
some additional information on the time evolution of the neck radius and the cavity. 
Finally, Appendix B discusses the dynamics of the minimal radius within the context of 
the model. 



2. Experimental and numerical results 

2.1. Experimental setup and procedure 

A sketch of the setup is seen in Fig. [T^. A disk of radius h^ is mounted on top of 
a thin rod (0 6 mm). This rod runs through a seal in the bottom of a large tank 
(500 mmxSOO mmxlOOO mm) and is connected at the lower end to a Thrusttube linear 
motor which is used to determine and control the velocity and acceleration of the disk. 
The position of the motor (and thus of the disk) along the vertical axis is measured 
with a spatial accuracy of 5 /im over a range of 1 m, the large acceleration of the motor 
(up to 30 5, with g the gravitational acceleration) makes it possible to perform impact 
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Figure 1 . (a) Schematic figure of the experimental setup with the linear motor controlling the 
rod and disk. The vertical rod runs through a seal in the bottom of the glass water tank and 
is pulled down by the linear motor in order to impact the disk on the water surface. Here, ho 
is the radius and V the controlled and constant velocity of the disk, (b) The formation of an 
upward and downward jet after the cavity has closed. In the present study we will focus on the 
cavity dynamics until pinch-off, just before the jet formation. 



experiments with constant velocities up to 5 m/s. 

The effect of the small diameter of the rod on the global flow and dynamics of the 
cavity is assumed to be negligible. As the minimum radius for the disk used in the 
experiments is 10 mm, the ratio of the cross-sectional area of the rod and the surface 
of the disk is always smaller than 9%. Since the rod is mounted in the center of the 
disk, where stagnation would normally occur, the influence on the radially outward flow 
below the disk is presumably small. 

Using the flat plate approximation, we can also estimate the direct contribution of 
the boundary layer of the rod to the axial flow. The boundary layer thickness d for a 
flat plate is given by Blasius' solution 6 w Syi^At, where v is the kinematic viscosity of 
water and At the time the boundary layer has to develop. We will equate the time At 
to the duration of the experiment, namely to the time interval starting from the impact 
of the disk until the collapse of the void, which in our experiments is found to scale as 
At sa 2.5-\//io/g, as will be discussed in detail in subsection 14. II For the largest disk size 
in the experiment ho = 40 mm, this result predicts a maximum boundary layer thickness 
of 1.8 mm. Under most experimental conditions of this study it is considerably thinner. 

In our experiments we pull the disk down with a constant velocity V. Making this 
main control parameter dimensionless, we obtain the Froude number Fr — V^ /{gho). 
The liquid properties are expressed in terms of the Reynolds number Re = Vho/v and 
the Weber number We = pV^ho/a, where a denotes the surface tension and p the 
fluid density. Since the Reynolds number and the Weber number are considerable on 
the large scales of Fig. [l] the viscosity and the surface tension do not seem to play a 
role. To be more precise, in our experiment the Reynolds number ranges between 500 
and 1.6 • 10^ and the Weber number ranges between 34 and 8.8 • 10'^. Note however 
that under only slightly different conditions, namely replacing the disk b y a cylinder 



submerged in water to avoid the splash, capillary waves do play a role (see iGekle et i 
(|2008 )). For the impact of a disk we find the only important dimensionless parameter 
to be the Froude number, i.e., the ratio of kinetic to gravitational energy, which 
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ranges from 0.6 to 200 in our experiments. It is convenient to use the amount of time 
r remaining until cavity collapse which is given by r = t—tcoii with tcoii the collapse time. 



2.2. Numerical method 



ine n u merical calculations are p e riormed usmg a bounaary in tegral metnoa (|rrosperetti 
(|2002h : IPower fc Wrobell (| 19951 ): lOguz fc Prosperettil (|l993l )) based on potential flow. 



This assumption excludes viscous effects and vorticity, which due to the short duration 
of the phenomenon and the high Reynolds number seems reasonable. 
Our code uses an axisymmetric geometry thus reducing the surface integrals to one- 
dimensional line integrals. For the time-stepping an iterative Crank-Nicholson scheme 
is employed. The size of each time step is calculated as t^tej} = f • iiiin(t„o£;e) with 
tnode = ^s/vnodei whcrc As is the distaucc to the neighboring node and w„ode the local 
velocity. With the safety factor / chosen to be 5% this procedure reliably prevents 
collisions of two nodes which would lead to serious disturbances in the numerical scheme. 
The number of nodes is variable in time, with the node density at any particular point 
on the surface being a function of the local curvature. This procedure guarantees that in 
regions with large curvatures, especially around the pinch-off point, the node density is 
always high enough to resolve the local details of the surface shape. At the same time, 
no computation power is wasted on an exceedingly high node density in flat regions 
towards infinity (which in our simulations is chosen to be 100 disk radii away from the 
central axis). To avoid numerical disturbances, we employ a regridding scheme in which 
at every second time-step the surface nodes are completely redistributed placing the 
new nodes exactly half-way between the old nodes. 

A particularly sensitive issue is the modeling of the crown splash created when the 
disk impacts the water surface. After first shooting upwards in a ring shape, the splash 
quickly breaks up into a large number of drops (which are ring-shaped due to the 
imposed axial symmetry). These drops do not further influence the cavity behavior 
and therefore need not be accounted for in our numerical code. In most simulations 
presented in this work, the crown-splash evolves normally until drop pinch-off. As this 
happens, the surface is reconnected at the pinch-off location and the drop is discarded. 



2.3. Interface 

The series of events typical for the experimental range of 1 < Fr < 100 is seen in the 
snapshots of Figs. [2^, b, and c. Upon impact a splash, an outward moving crown of 
water, is formed. A void is created which collapses due to the hydrostatic pressure and 
a singularity arises when the collapsing walls of the void collide with each other. Two 
jets emerge in this experiment: One upwards straight into the air, and one downwards 
into the entrained air bubble (see Fig.[T]3). 

In each of Figs. [2^, b, and c the experimental sequence is overlaid with the results of 
our boundary integral simulation. For Fr=0.85 and Fr=3.4 (Figs. [2^ and b), the cavity 
dynamics is found to be captured extremely well by the numerical result, represented 
by the overlaid lines. Note that this is a one-to-one comparison between simulation and 
experiment, without any rescaling in time or space. Due to the axisymmetry of the 
code it is not possible to capture the full details of the splash and since our focus is on 
the cavity dynamics we chose to simply take out any droplets which are released from 
the splash. Surface tension however still expresses itself in small capillary waves in the 
region of the splash. These waves are most notable in Fig. [2^. As was mentioned before, 
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Figure 2. Snapshots of the formation and collapse of a surface void in the plunger experiment: 
A linear motor pulls down a disk of radius ho — 30 mm through the water surface at a constant 
velocity (a) V = 0.5 m/s (Fr = 0.85); (b) V = 1.0 m/s (Fr = 3.4); (c) V = 2.0 m/s (Fr = 
13.6); and (d) Fr = 200. The collapse of the void is imaged at 1000 frames per second. The lines 
(overlay) are the void profiles obtained from boundary integral simulations. Without the use of 
any free parameter, neither in time nor in space, an excellent agreement between the simulation 
and experiment is found in (a) and (b). Due to a (mild) surface seal there is a discrepancy 
between the simulations and the experiment in sequence (c), both in the top region near the 
splash and in the pinch-off region. The region of the dashed box is shown enlarged in Fig. |3l 
The experimental data for (d) (not shown) is severely dominated by a strong surface seal, which 
is an air effect and not the focus of the present study. 
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Figure 3. Enlargement of the region around the pinch off point at r = 1 ms from the sequence 
Fig. [2}; (Fr — 13.6). A significant discrepancy can be seen for the depth of the pinch-off between 
the boundary integral simulation (white line) and the experimental recording. The origin of this 
discrepancy is the airflow in the cavity as will be elaborated elsewhere. 
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Figure 4. Snapshots of the surface seal which occurs for a disk of radius ho = 30 mm impacting 
the water surface at a constant velocity of V^ = 2.0 m/s (Fr = 13.6), i.e., under the same 
conditions as Fig. [5J;. 



similar capillary waves (but from a different origin) are found to have a si g nifica nt 



influence on the closure of the cavity for a submerging cylinder ( Gekle et all ( 20081 )). 



For the impacting disk discussed in this paper however they do not affect the closure. 



The results for Fr=0.85 in Fig. ^ji. illustrate the effect of the relative importance of 
gravity. In the last frame of Fig. [2^ it can be seen that the cavity is less symmetric 
in axial direction around the closure point compared to the experiment performed at 
Fr=3.4 shown in Fig. [2)3. In the third sequence at Fr =13.6 (Fig. [Sfc), which goes beyond 



the experimental Froude number range described in lBergmann et ali { 200(f ). significant 
deviations between the experiment and the numerical cavity shape are found, most 
notably in the enlargement of Fig. [3] at the depth of the cavity closure. The closure of the 
cavity is found to be somewhat deeper in the numerics as compared to the experiments. 
This deviation can be attributed to a secondary effect due to the su rrounding air , called 
the surface seal (see Fig. |4|). This phenome non was first reported bv W prthingtonI ( 19081 ) 
and later investigated in more detail by iGilbarg fc Anderso n (1948). Note that the 
impact experiment of Fig. [5] is performed under the same conditions as Fig. ^. The 
surface seal is the entrainment of the initially outward moving splash by the air rushing 
into the expanding cavity. If the airflow is strong enough, the splash will close on the 
axis of symmetry and completely seal off the top of the cavity above the height of the 
undisturbed water surface. 



The surface seal is found to become more pronounced at higher impact velocity, where 
the surface seal occurs earlier and more liquid is involved in this closure. Accordingly, 
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there is also a larger influence on the shape of the cavity at higher impact velocity. Since 
this article aims to treat the purely pressure driven collapse of the cavity, without the 
contributions of the surrounding air, our experimental range is limited by the occurrence 
of the surface seal. In the simulations we therefore intentionally do not include the 
air. This explains the discrepancy of Fig. [2}; (enlarged in Fig. [3]), since contrary to the 
experiments, no surface seal occurs in the numerics due to the absence of air. In Fig. [2Ji 
we go far beyond the experimentally available range by performing simulations at a 
Froude number of 200. 

It is instructive to compare t he present boundary integral simulation results with 
those reported by iGaudetl ( 1998 ). who reported a bulging contraction of the cavity at 



the surface level. He found this contraction to close for Fr > 200 and interpreted it as 
a surface seal in the absence of air. We found no evidence for such a surface seal in our 
simulations, even fo r cons iderably larger Froude numbers, and surmise that the effect 
reported bv iGauded ( 19981 ) may be connected to using an insufficient number of nodes 



in the splash region caused by the limited amount of computational power available at 
that time. 



2.4. Flow field 

In the previous subsection we found the experimental shape of the impact cavity to 
be well described by our boundary integral simulations if no surface seal occurs. The 
question we will address in this subsection is whether the simulations also give an 
accurate description of the surrounding flow field. To this end we will measure the 
velocity field around the transient cavity through high speed particle image velocimetry 
(PIV). These experiments are crucial to check the validity of the boundary integral 
simulations, as the presence of a solid boundary, namely the impacting object, will induce 
vorticity in the flow. We will compare the experimental flow field to the boundary in- 
tegral results and finally investigate the radial flow at the depth of closure in more detail. 

To perform the PIV measurements, the fluid is seeded with small DANTEC Dynamics 
polyamid tracer particles of radius 25 /im and density 1030 kg/m"^ which follow the 
flow. A laser sheet shines from the side through the fluid, creating an illuminated 
plane through the symmetry axis of the cavity. The light scattered by the particles 
is captured by a high speed camera at a frame rate of 6000 frames per second and 
a resolution of 1024x512 pixels. The series of recorded images is then correlated by 
multipass algorithms, using DaVis PIV software by LaVision, in order to determine 
the flow fleld in a plane in the liquid. The correlation was performed in two passes at 
sub-pixel accuracy, using 64x64 pixels and 32x32 pixels interrogation windows. The 
windows overlap by 50%, resulting in one velocity vector every 16x16 pixels. 

In order to obtain high resolution PIV measurements of the flow around the cavity, 
we made use of the reproducibility of the experiment. The left side of each of the images 
of Fig. [5] shows the flow around the expanding void by combining the results of four 
separate PIV measurements at different depths. In this fashion PIV experiments were 
performed for a field of view of 96 mm x 56 mm at a spatial resolution of 0.9 mm (In 
Fig. [5] only 0.7% of the measured vector field is shown). This high resolution makes it 
possible to simultaneously compare the global fiow, as well as the smaller fiow structures 
at the pinch-off depth and the disk's edge. 

The right side of each image of Fig. [5] shows the numerically obtained cavity profile 
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Figure 5. Comparison of tlie flow fleld obtained from experiments and boundary integral sim- 
ulations for a disk of radius hdisk ~ 20 mm which impacts with a velocity oi V — 1.0 m/s (Fr 
= 5.1). The flgures show the flow fleld at r = (a) 49 ms; (b) 25 ms; (c) 7 ms; and (d) 1 ms. The 
left side of each image shows the flow field (overlaid vectors) obtained from the experiment by 
four high-speed PIV recordings at 6000 frames per second. The four separate recordings were 
taken at different depths and combined to give the flow fleld at high resolution. The recordings 
on the left side of each image also illustrate the degree of reproducibility of the experiment, as 
the match between the four PIV recordings at different depths obtained from four repetitions 
of the same experiment (see main text) is perfect. For clarity only 0.7% of the measured vectors 
is shown. The right side of each of the images shows the void profile and the corresponding flow 
field (overlaid vectors) obtained from the boundary integral calculations. 



and flow field. At first siglit tliere appears to be a good agreement, but one would like 
to obtain a more quantitative comparison between experiment and simulation. This is 
provided in Fig. [SI which shows contour plots of the axial flow component (Fig. [S^-d) 
and the radial flow component (Fig. ^-f) obtained from the PIV measurements (at 
the left side of each image) and boundary integral simulations (at the the right side 
of each image). From this figure it is clear that the magnitude as well as the topology 
of the flow are in excellent agreement. Figures [5] and [6] are a one-to-one comparison 
between simulation and experiment, and we stress once more, without the use of any 
free parameter. 

In addition to the above, the experimental pictures of Fig. [S] reveal that our initial 
assumption to neglect the influence of the rod on the flow (see Section [2]) is correct. 
The rod itself is clearly visible in the experimental snapshot of Fig. [5^ and the PIV 
software has correctly detected its downward movement, as can be seen in Figs. [5^-b. 
From the same figures we also conclude that outside a thin region around the rod the 
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Figure 6. The axial (left figures) and radial (right figures) components of the fiow field, repre- 
senting the same velocity data as in Fig. [5] and taken at the same times t — 49 ms (a,e), 25 ms 
(b,f), 7 ms (c,g), and 1 ms (d,h). The four images to the left (a-d) compare a contour plot of the 
axial flow component from the experiment (left side of each image) with that of the numerics 
(right side of each image). The four images on the right (e-h) show a similar comparison in a 
contour plot of the radial flow component from the experiment (left side of each image) and 
numerics (right side of each image). Apart from the region where the rod is pulling down the 
disk in the experiment, which is absent in the simulation, both components of the flow field 
show excellent agreement between the experiments and numerical calculations. Again this is a 
one-to-one comparison of experimental and numerical data without any free parameter, neither 
in velocity, space, nor time. 
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Figure 7. Radial (a) and axial (b) component of the flow at the depth of closure Zcoii for a 
disk of radius hdisk ~ 10 mm which impacts at a velocity of \^ = 1.0 m/s (Fr = 10.2). The 
symbols show the result obtained from the PIV measurements at different times, () : t — 30.5 
ms; D : r = 20.5 ms; V : t = 10.5 ms; and o : t = 0.5 ms, and the solid red lines are the 
corresponding numerical results from the BI simulation. Note that these lines end on the cavity 
surface (solid red dots). The PIV data is an average of 6 subsequent measurements obtained 
from the high speed PIV recordings at 6000 frames per second. In consequence, the Vr and Vz 
velocity components shown here are the average over one millisecond. 



flow remains unchanged. Most importantly, the outward flow at the edges of the disk in 
Fig. [6^, which is responsible for the expansion of the void, is unaffected by the presence 
of the rod. This can be understood from the fact that below the disk the radial flow 
component decays quickly towards the center of the disk whereas the vertical component 
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in the center is equal to the disk speed. As the fluid in the central region hardly moves 
with respect to the disk, the presence of the rod has a negligible eS'ect on the flow field. 
We conclude that the boundary layer of the rod makes no contribution whatsoever to 
the flow field around the cavity. 

In Fig. [7] the experimental flow at the closure depth is compared with numerics up 
to 7 disk radii in radial direction in order to obtain a more quantitative measure of the 
magnitude of the deviations between the numerical and experimental flow field. We find 
this deviation to be typically of the order of 0.01 m/s, but it can be slightly larger if the 
flow velocity is small. This larger inaccuracy at low flow velocities is generic to the PIV 
method and can most clearly be seen to occur for r — 10.5 ms in Fig. [TJs. Overall, a 
very good agreement is found between the far field flow in the numerics and experiments. 

Both in the experiment and simulation we observe that during the expansion of the 
void the magnitude of the outward radial flow falls off with the distance to the symmetry 
axis (Fig. [7^). However, once the cavity starts to collapse inward there will be a region 
around the cavity where the (radial) direction of the flow is reversed and there will be an 
axisymmetric curved plane (manifold) where the radial flow component vanishes. Here 
this happens between r = 30.5 ms and r = 20.5 ms (cf. Fig. [7^). In subsection 13.21 we 
will discuss in detail how this reversal of the radial flow expresses itself in the radial 
dynamics of the cavity. 



3. Modeling the cavity dynamics 

In this section we will first derive a simple analytical model for the radial dynamics of 
the transient cavity. Secondly, we will investigate the surrounding flow, which enters the 
model through two of the free parameters and causes an asymmetry of the collapse. In 
the last part of this section we compare the model to the radial dynamics of the cavity 
observed in experiment and simulation. 

3.1. A model for the radial cavity dynamics 

The full analytical modeling of a cylindrical symmetric collapse of the transient cavity 
presents the difficulty of a coupling between the free surface and the flow surrounding the 
cavity. To tackle this difficulty we propose the convenient simplification of dividing the 
problem up into a set of quasi two-dimensional problems. If the axial component of the 
flow is small compared to the horizontal flow components, we can approximate the flow 
as to be confined to the horizontal plane. In this way an equation for the collapse of a 
two-dimensional bubble will suffice to describe the cavity dynamics at an arbitrary depth. 

To derive such an e quat ion we wil l clos ely follow a derivation given in 
Qguz fc Prosperettil (|l993l ) and lLohse et al\ (|2004l ). The argument starts by writing the 



Euler equation in cylindrical coordinates, thereby neglecting the vertical flow component 
and its derivatives. This means that we assume the flow to be quasi two-dimensional at 
any depth along the cavity. The azimuthal components can be ignored due to the axial 
symmetry, leaving the following equation 

dvr dvr I dp 

-^+'^r^^ = — ^, (3.1) 

at or p or 

where p denotes the density of the liquid. Under the above assumption of negligible Vz 
and thus dvz/dz, the continuity equation and the boundary conditions on the surface of 
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Figure 8. Schematic representation of the three stages in the radial dynamics of the cavity at a 
constant depth z. The first two stages (corresponding to the time interval Atexpa + Atctra) are 
governed by a forcing of the flow by hydrostatic pressure. We can distinguish an expansion and 
a contraction phase for which the model parameter hoo differs considerably. In the third stage 
(corresponding to Atcoii) the collapsing void accelerates towards the singularity (pinch-off) in 
which inertia takes over as the driving factor. 



the void lead to a second equation 



rvrir,t) = h{t)h{t). 



(3.2) 



Here, h{t) is the radius of the cavity and its derivative h{t) the velocity of the cavity wall 
at any depth z below the surface. Substituting Eq. (|3.2p into Eq. (|3.ip gives 




= 0. 



(3.3) 



We can integrate this equation over r from the cavity wall ft, to a reference point hoo, 
where the flow is taken to be quiescent. This integration yields a Rayleigh-like equation 
for the void radius at a fixed depth z, 



d{hh) 
dt 



log- 



1 



gz. 



(3.4) 



Here, we have used the fact that the pressure (poo) driving the collapse of the cavity is 
provided by the hydrostatic pressure pgz, where z denotes the depth below the fluid 
surface. Close to the collapse, the quantity hoo can be interpreted as the length scale 
related to the matching of an inner and outer flow region. In the (inner) region near 
the neck the induced flow looks like a collapsing cylinder as described by Eq. p.4p . 
whereas in the (outer) region far from the neck, the flow resembles that of a dipole 
(plus its image in the free surface). A complete description of the flow would require 
the matching of these two regions, where hoo would be determined in the process as the 
cross-over length scale, hoo would thus be expected to be of the order of a typical length 
scale of the process, such as the distance of the cavity surface to the plane where the 
radial flow vanishes (see Fig. [7^). Therefore, strictly speaking, hoo is a function of the 
Froude number and time. In the model presented below we follow a different, simplified 
route and set hoo to a constant value (a time averaged value of its dynamics). 



We will now use Eq. 



to analyze the radial dynamics from the initial impact of the 
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disk to to the time of closure of the cavity tcoii at arbitrary depth z. In order to obtain an 
analytical approximate solution, we decompose the cavity dynamics into three different 
stages, depicted schematically in Figure[Hl In this figure the time intervals corresponding 
to the different stages are denoted as Atg^pa, Atctra, and Atcoii respectively. In the first 
two stages, during At^xpa and Atctra, the dynamics is dominated by the hydrostatic 
pressure forcing. In these stages we observe that the water is first pushed aside by the 
passing disk, creating an expanding void. At the maximum radius h^nax, the expansion 
is halted and the void starts to contract, h^ax is typically of the order of Hq, e.g., for 
Fr = 3.4 and Fr = 200 we find respectively hmax ~ 1.3/io and hmax ~ 2.4/io- Since 
hifmax) — 0, we can assume that h{t) is small during this expansion and contraction and 
we can neglect the second term in Eq. p.4|) leading to 

djhh) 
dt 



log -— =gz. (3.5) 



Since log{h/hoo) varies very slowly in the first regimes, we equate log{h/hoo) ~ 
log{h(tmax)/hoo) and we solve Eq. ([XS)) using hit^ax) = h^ax and h{tmax) = 0, leading 
to a parabolic approximation for /i^, 

h'^{z,t) = hl^^^- —{t-trnaxf-, (3.6) 

P 

with /3 = ~-\og{h^nax/hoo)- The above equation holds for both the expansion stage, the 
time it takes for the void to grow from h^ to hmax, and the contraction stage, the time 
it takes to shrink back to Hq. 

In the third stage, during Atcoii, the collapsing void accelerates towards the singularity 
and inertia takes over as the dominant factor driving the dynamics of the cavity. This 
stage can be described using a different approximation to Eq. p.4p . Near the collapse, 
h approaches zero, /loo is typically very large and thus the logarithm diverges. The only 
way Eq. (|3.4p can remain valid is when the prefactor of the logarithm vanishes. This 
means that 

djhh) _l d^h') _ 

Integration gives the power law of the two-dimensional Rayleigh collapse (cf. 
iBergmann et aU (|2006i) ) 

h{t) = ^/C{t,oU-tf'^. (3.8) 



In subsection 13.31 the integration constant yC will be determined from continuity of h 
and h. 

3.2. The influence of the flow on hoc 

As an intermezzo in the exposition of the model we now turn to an important point, 
namely that there is a significantly different quality to the flow in the expansion and 
the contraction stage. This difference already became clear in our discussion of Figure [7] 
where we found that in the expansion phase the outward radial flow simply decays with 
the radial distance, whereas in the contraction phase the radial flow becomes zero at 
some finite distance at which it changes direction. This is due to the fact that the fluid 
flows outward until the cavity reaches its maximum radius h^ax , from where it will start 
to move inward, creating a reversed-flow region around the cavity wall which grows in 
time. Although in both stages hydrostatic pressure is the dominant factor driving the 
dynamics of the cavity, there is this dissimilarity in the surrounding flow which needs 
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Figure 9. (a) The radial velocity component Vr at the moment and depth of closure multiplied 
by the radial distance r as a function of the same radial distance for different Fronde numbers. 
The inset shows the original radial velocity data. The distance where the radial velocity is zero 
and the flow is stagnant in radial direction is indicated by the blue stars. The flow can be seen to 
resemble a radial sink flow more with increasing Fronde number as the minimum of VrV decreases 
with increasing Froude number. 

(b) The distance of the point where the radial flow reverses sign to the symmetry axis, determined 
at the depth of closure as a function of the normalized time remaining till closure r = tcoii — t. 
Below the curves the flow is directed inwards, above them it is directed outwards. The radial 
distances of the flow reversal point at r = in this figure correspond to the blue stars of figure 
(a). The distance of the point of flow reversal is related to the length scales ft^'^" and /i^*"", 
which are therefore expected to behave similarly in time. 

to be incorporated into the model. To investigate the dynamics of this dissimilarity 
in detail, we turn to the simulations from which we can obtain the flow field with an 
arbitrarily fine resolution. 

Figure [9K shows the radial flow component Vr multiplied by the radial distance r to 
the axis of symmetry at the depth and moment of pinch-off. Since the flow at the neck 
resembles locally a two dimensional sink, whose strength falls off with l/r, multiplying 
Vr with r eliminates this geometrical contribution to the flow. For the lower Froude 
numbers of 3.4 and 5.1 the radial flow component reverses direction at closure depth 
and time at some distance r (blue stars). At the higher Froude numbers (10.2 and 200) 
no such point is observed within the numerical domain, which extends to 100 disk radii 
in the radial direction. This docs not mean that such a flow reversal point is absent 
during the complete time of the collapse, as can be seen in Fig. [9Jd where we plot the 
location of the flow reversal point as a function of time: The radial flow reversal point 
comes into existence at the wall of the void at the moment that the expanded cavity 
starts to collapse and the flow direction is reversed inward. From then onwards, this 
point travels away from the axis of symmetry as the collapse is approached (r ^ 0). In 
the same figure we also observe that for a higher Froude number the radial flow reversal 
point travels outward much faster during the cavity collapse as compared to the low Fr 
case. 



The reversed-flow region can be characterized by a stagnation point (or saddle point) 
in the (r,z)-plane, corresponding to a circle in three dimensions, where both the axial 
and the radial velocity components change sign. In Fig. (TU] the path of this stagnation 
point is shown for Fr — 3.4 (a) and Fr = 10.2 (b). For both simulations the stagnation 
point not only moves away from the axis of symmetry as the pinch-off is approached, 
but it is also seen to move down in the axial direction and at some point even to cross 
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Figure 10. The open circles show the path of the stagnation point of the flow for Froude 3.4 
(a) and 10.2 (b). For (a) the first observation of the stagnation point is made at r = 61 ms until 
1 ms before closure at intervals of 3 ms. In (b) the flrst observation is made at r = 35 ms until 
7 ms before closure at intervals of 1 ms. For clarity the time till closure is indicated only for 
every second observation. Further more, the void profiles at the time of the first (blue) and last 
(red) observation are shown. The depth of closure is indicated by the dotted line. 



the depth of closure. A similar path of the stagnation point is observed for all the 
simulations of Fig. [9] and only at one instant during the collapse of the cavity does the 
radial flow reversal point at closure depth truly coincide with the stagnation point. 



The above leads us to three observations which are relevant for our model of the 
cavity collapse: (i) Since radial flow reversal at closure depth occurs when the cavity 
starts to collapse, the topology of the flow differs between the expansion and contraction 
stage. Since hoc is the radial distance at which the flow can be assumed to be quiescent 
(v = and p = poo = pgz) it is related to the structure of the surrounding flow, and 
it is therefore justified to assume different values of ft-oo in the respective stages. We 
will take hoc = ^^^° in the expansion stage and hoo = h'^'^ in the contraction stage of 
the model. Just like hoo in subsection 13.11 h^P"" and /ij^'"" are set to a constant value, 
representing the time averaged behavior of hoc in each respective stage, (ii) As the 
distance of the radial ffow reversal during the contraction moves away faster at higher 
Froude number, presumably a higher val ue for h'^"' needs to be taken for larger Froude 
numbers, (iii) In iBergmann et al\ (|2006f ) we found that there are two scaling regimes 
for the neck radius, a first regime where the neck radius scales as a pure power law of 
time (as in Eg. 13. 8p . and a second regime, where a logarithmic correction of time has to 
be taken into account. The crossover between both regimes is given by the length scale 
/^maa/^ocT"- As we find from Fig. [5]d, for all Froude numbers the distance of the radial 
flow reversal increases as the pinch-off is approached. Although in theory we assume 
h'^'^ to be constant, in reality h'^"' thus increases as the pinch-off is approached. This 
means the cross-over length scale h'^^^/h'^"' decreases with time. 

Therefore the time needed for the collapsing neck to decrease to h'^^^/h'^"' will be 
longer as compared to the assumption of a constant (initial) value for h'^'^ and may 
even never reach this second regime. The effect is stronger for increasing Froude number, 
since the radial flow reversal point at closure depth moves away faster and further at 
higher Froude number. 
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3.3. The free parameters of the model 

In this subsection we continue our derivation of a simplified model for the radial cavity 
dynamics started in subsection l3.1l As argued in the previous subsection it is justified to 
assume different (constant) values for h^o during the expansion and contraction stage of 
the void. We therefore introduce different values for (3 in Eq. (|3.6p . depending on whether 
we are in the expansion or in the contraction stage 

Q J Pexpa ^ ^^^y'^maxy^ ) / '^OQ,expa) ^ ^ ^max /o q\ 

^ Pctra ^ ^^hy'^maxy^ ) / '^OQjCtra) ^ ^ ^max 

Note that with this definition Pexpa and Pctra are positive quantities as for both holds 
hoc,expa,haa,ctra 3> h^nax- Secondly, the fact that /3 depends only logarithmically on 
/loo furthermore justifies approximating the time-dependent quantity hao(t) by its time- 
average hoc- 

Now, to determine hmax , or rather the time it will take to get there from the time the 
disk passes at i = treach, we need the radial velocity of the initial expansi on at t = treach 



(see F ig. [8l . A reasonable assumption (and similar to the proposition of lDuclaux et 



(120071)) is that the disk displaces water from underneath itself to the sides at a velocity 



directly proportional to its downward velocity. Therefore, we have 

h{treach) = OlexpaV . (3.10) 

For the velocity at the end of the contraction phase at t cross = treach + Ateipa + ^tctra, 
we write in a similar fashion 

h{tcross) = -OictraV . (3.11) 

Clearly, both Uexpa and Udra are again positive quantities. 

The analytical model for the radial cavity dynamics given by Eq. (|3.6p and Eq. (|3.8p 
thus has four unknown parameters Uexpai Pexpa i oidra, and Pctra- The value of C in 
Eq. (j3.8p follows from the fact that the trajectory and its derivative in Fig. [8] should be 
continuous. We assume the collapse regime to start at the end of the contraction phase, 
where we have h{tcross) — ho and h{tcross) = —o^ctraV. From these conditions, the value 
of C is readily obtained, 

C^2hoactraV. (3.12) 

However, for given aexpaj Pexpa, and Pctra, the constant adra is also uniquely determined 
by the continuity of the trajectory and its derivative at h{tjnax) — hmax j which gives 

^ctra — ^expa\J Pexpa / Hctra -i yo.LO) 

and leaves only aexpa, Pexpa, and Pctra to be determined. 

Summarizing, the time evolution of the cavity at depth z is described by the following 
three equations 

h{z,t) = Jhl^^.-^ - J^it - traaxf- iOT treach < t < t^ax , (3.14) 

fl[Z,t) ^= \/ il-rnax Batra ^ ''max) tOr tmax < t S tcross , (o.ioj 

h{z,t) = V^haactraV^/tcoU - t iov tcross < t < tcoll , (3.16) 

where the times treach, tmax, tcross, and tcoii are readily related to the impact time t — 
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Figure 11. (a) The time evolution of the radius of the cavity at closure depth for different 
Froude numbers. The solid blue lines represent the simulation results and the dashed red lines 
correspond to a least square fit of the approximation given by Eqs. (|3.14|l - (|3.16p . (b) The same 
numerical ti me evolution d ata as in (a) (solid blue lines) are now approximated by the model 
proposed bv lDuclaux et al\ ([2007), which consists of Eqs. (|3.14|l - (l3.15p with /Sexpa and Pctra set 
to 1, again from a least squares fit (dashed red lines). In both (a) and (b) time has been rescaled 
by a factor Vho/gzcott in order to show the results for the two Froude numbers in a single plot. 



(which will be done explicitly in section [4]) and hmax is given by 



hmax (z) = hoJl + al /3, 



parexpa 






(3.17) 



as can be easily derived, e.g., from Eq. 13.141 together with its boundary conditions 

h{z,treach) = ^0 and h{z,treach) = CtexpaV . 

3.4. Validation of the model 

We will now compare the dynamics of the radius of the void at closure depth with the 
theoretical prediction of Eq. (|3.6p and Eq. (|3.8p to validate the model and quantify the 
influence of the flow reversal on Pexpa and Pctra- 



The parameter adra is eliminated by the relation Eq. (|3.13l) . leaving three parameters 
to match Eqs. (|3.14p - (|3.16p to the radial dynamics from the boundary integral simu- 
lations. At first sight, one could assume that the initial outward velocity aexpaV could 
be easily obtained from simulation or experiment, since it is observed as the angle at 
which the free surface leaves the disk. However, when this angle is investigated in closer 
detail, it is found to strongly depend on the distance from the disk's edge over which 
it is measured. In the numerics, close enough to the disk's edge, the free surface even 
becomes nearly parallel with the disk. This means that although aexpoV is useful as a 
(theoretical) concept, it is not directly measurable and should therefore be determined 
by a fitting routine. 

The three parameters aexpa, Pexpa, and Pctra are determined by a least square fit 
to the radial dynamics of the cavity at closure depth obtained from the simulations. 
Fig. [TTb shows the comparison between these fits of Eqs. p.l4p - p.l6p (red dashed line) 
and the simulations (blue solid line) at two different Froude numbers of 3.4 and 200. The 
approximation is found to be in excellent agreement throughout the collapse, faithfully 
reproducing the maximum expansion of the cavity and the complete time of collapse. 

In Figure [T^ we find the parameters otexpa, Pexpai and Pdra as a function of the 
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Figure 12. (a) The parameters a^^pa, Pexpa, Q-ctra, and fictra obtained from the fitting rou- 
tine used in Fig. Illb .. All are found to weakly depend on the Froude number (note the loga- 
rithmic scale for Fr). (b) The quantities 'J OLexpaii^x-pa + actraf3ctra and ala:paPexpa calculated 
from (a). Whereas y/aexpaPexpa + ctctraPctra is found to weakly depend on the Froude number, 
a 



Pe^ 



expaHexpa 



has a nearly constant value of 0.40. (horizontal dotted line). In comparison with 



\/oLexpaiiexpa ■\- OctraPctra, the horizoutal dashcd line indicates the value Ci = 1.10 obtained 
from a best fit to the closure depth data of Fig. 1131 



Froude number, determined by repeating the fitting routine described above for many 
impact velocities. All are found to weakly depend on the Froude number (note that 
a logarithmic scale has been used for Fr). For completeness we also plot the derived 
quantity actra, calculated from Eg. 13.131 



If in Eqs. p.l4p ~ p.l5p the constants f3expa and Pctra are set to 1 an d therefore by 
we arrive at the cavity dynamics proposed by iDuclaux et al. 



O^ctra 



^expa: 



(|2007f ) for impacting spheres and cylinders. These dynamics are shown in Fig. Illb with 
the only free parameter aexpa also determined by a least square fit to the data. This 
approximation is seen to qualitatively reproduce the trend for the maximum expansion 
an d collapse time , but fa ils to capture the exact values. It is also conceptually different, 
as 



Duclaux et 



( 20071 ) propose the cavity dynamics to be symmetric around the 
maximum expansion, while our model captures the asymmetry around this point in 
time that is also found in experiment and simulation. Our solution Eqs. (|3.14p - (|3.15p is 
explicitly not symmetric, since it allows for different values of /3 at t < t^ax and t > tmax- 



To conclude this section we return to the first two observations we made at the end 
of subsection 13.21 on the motion of the stagnation point and the plausible consequences 
for hoo- (i) The flow reversal which occurs when the cavity starts to collapse indeed 
introduces an asymmetry in the behavior around the maximum expansion. This is 
clearly observed in the radial dynamics of Fig. [Til especially for Fr = 3.4. (ii) As the 
distance of the radial flow reversal during the collapse moves away faster at higher 
Froude number (see Fig. ^p) , we indeed have to introduce a larger ft.^'"" (corresponding 
to a larger Pctra) for higher Froude number in the fit of Fig.fTTk to account for this effect. 
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4. Characteristics of the Transient Cavity 

Now that we derived a simplified model for the radial dynamics of the cavity, we will 
use it, together with the simulations and experiments, to investigate the following key 
characteristics of the transient cavity: (i) the depth of the pinch-ofF and the depth of the 
disk at the moment of pinch-ofF (subsection 14. ip , and (ii) the amount of air entrained 
through the cavity collapse (subsect ion 14. 2p . Besides this we elaborate on our previous 
findings from lBergmann et all ( 2006[) in the Appendix [Xj There we revisit the dynamics 



of the cavity at closure depth (Appendix lA.ip and the cavity shape around the minimal 
radius (Appendix IA.2[) . Finally, in Appendix [B] wc discuss the time evolution of the 
vertical position of the minimal cavity radius and place it within the context of the 
model. 

4.1. Closure depth 
Following [Clasheen fc McMahonI (|l996l) . iGaude^ (|l998l ). and iDuclaux et all (|2007l ) we 



will characterize the shape of the cavity at pinch-off by the depth of closure Zcoii, i-e., 
the depth at which the pinch-off takes place. To capture more information on the full 
shape of the void, we will also investigate how Zcoii relates to the total depth of the 
cavity zj^isk{tcoii) = Zdisk,coU at the time of collapse (or closure) (see the inset of Fig. [T^ 



A comprehensive argument for the scaling of Zcpu can be obtained by following a 



similar procedure to the one outlined in iLohse et al\ (|2004r ) for the determination of the 



closure depth after the impact of a steel ball on soft sand. The difference is that whereas 
in sand one can assume that due to the compressibility of the material there is hardly 
any outwards motion of the sand, here we are dealing with an incompressible fluid and 
the outward expansion of the cavity needs to be taken into account. 

The time interval between impact of the disk and collapse of the cavity At — tcoii — to 
at any depth z consists of two main parts: First, the disk needs an amount of time Atreach 
to reach the depth z. Second, just after the disk passes there is the time Atyoid it takes 
for the void to form, expand, and collapse 

At = Atreach + Atyozd ■ (4.1) 

The first term equals Atreach = z/V since the velocity of the disk is constant in the 
experiment and simulation. In Section [3] Atyoid was decomposed into three stages as is 
schematically depicted in Fig. [51 The collapse time can thus be written as. 

At = Atreach + Atexpa + Atetra + Ateoll ■ (4-2) 

" V ' 

To estimate these last three timescales at arbitrary depth 2, we turn to our model for 
the cavity dynamics Eqs. (|3.14p - (|3.16p . 

If we combine the conditions Eq. (|3.10|) and Eq. (jS.lip with the time derivative of 
Eqs. (|3.14p and (|3.15p . we readily obtain 

A, - R ^ 

gz 

Atetra — CtetraPctra ■ (4-3) 

gz 

(4.4) 
Recall that aetra = aexpaiPexpa/ PetraV^^ ■ The radial collapse during Ateoii is in turn 
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described by the approximation of Eq. ()3.16p . Since hitcross) — ho we find for this time 
interval 

^tcoU — tcoll — tcross — 7; 77 • (4.5) 

^actra V 

Collecting all the above time intervals, within the model the total amount of time that 
passed from the impact of the disk until the collapse of the cavity at depth z is given by 

_£ , / « _L n ^/iQ^ , 1 ho 

T^ ' [P^expaPexpa 1 O^ctraPctra) ~r ■ l^-*-*/ 

V gz 2actra V 

Now, to find the closure depth Zcoiu we need to determine at what depth the collapse 
will occur first, which we can do by solving 

f^ = 0. (4.7) 

dz 

This gives 

r — y (^expaPexpa ~r (y-ctraPctra^^ • V ^^J 

tlo 

In addition, the total depth of the disk at the time of collapse, Zdisk{tcoii) = Zdisk,coiu 
can be obtained by inserting Eq. (|4.8p into Eq. (|4.6p to give Zdisk, coii = VAt, or 

"", ' ''° = 2y/aexpaPexpa + actraPctraFl^ + 7: . (4.9) 

When we compare these expressions with the experiments without a surface seal 
(blue symbols) and the numerical calculations (black circles) in Fig. [T5b we find a 
very good agreement with the prediction of Eq. (|4.8p . A fit to the data of Zcoii gives 
Zcou/hci — CiFr ' , with Ci — 1.10. The agreement of the experiments in which a 
surface seal occurs (red symbols) deteriorates for a fixed disk size with increasing Froude 
number, since the surface seal becomes more disruptive at higher impact velocities. 
In the same figure we find the experimental and numerical results for the total depth of 
the disk at closure Zdisk,coii- From the apparent power- law scaling it is clear that the 
constant l/(2actra) in Eq. (|4.9p has no significant contribution. The total depth of the 
void is found to scale as Zdisk,coii/hi:, = C2Fr^", with C2 = 2.49 close to the expected 

value of C2 = 2Ci = 2.2 that follows from Eqs. ((48)) and (|49|) . The fact that the closure 

1 /'? 
depth and the total depth have the same power-law scaling Fr ' indicates that the time 

from the initial impact of the disk to the time of closure of the cavity does not depend 

on the velocity of the impactiJ, since Ai = Zdisk, coii/V — C2\/ho/g- 



This is in agreement with the findings of iGlasheen k, McMahonI (19961), who ex- 



perimentally observed a similar scaling for the impact of a disk on a water su rface, 



although with a slightly different prefactor of C2 ~ 2.3. iDuclaux et al\ (|2007t ) also 
found the scaling of At (ex ^yhojg) for impacting spheres and furthermore reported 

Zdisk, coll/ho = '^Zcoii/hf) cx Fr ' in agreement with our observations. 

To investigate the data of Fig. [T3k more closely it is convenient to take the ratio of 
Zcoii/ Zdisk, coll (see Fig. 113b). According to Eq. (|4.8p and Eq. (|4.9p . this ratio should scale 



f Within the model, keeping the constant term in Eq. (|4.9p is equivalent to keeping the last 
term in Eq. (|4.6|) which would add a l/V-dependence to the closure time, vanishing for high 
Froude number. 
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as 

Zcou/zd^sk, coll = 2 + Fr"^/'^ w 2 (4.10) 

^O^ctrayJ O^expaPexpa i O^ctraPctra 

in the limit of large Froude number. In Fig. 113b the ratio of Zcoii/ Zdisk,coii hi the 
experiments without a surface seal (blue symbols) and the numerical calculations (black 
circles) are indeed close to the constant value of 2 (dashed black line), but at lower 
Froude number it decreases slightly contrary to the proposed scaling by the second term 
in Eq. ()4.10p . Although the second term of the ratio of Eq. ()4.10|) should become relevant 
when the Froude number is considerably small, this is not observed in Fig. 113b . This 
can be understood by noting that in the limit of small Froude number our assumption 
of non-interacting fluid layers from Eq. (|3.4p breaks down as gravity becomes more 
important and thus Eq. (|4.10|) is no longer valid. In Fig. [T3b it is again illustrated that 
the experiments with a surface seal (red symbols) deviate more and more from the 
simulations without air as the Froude number increases. 



The fits to the trajectories discussed in subsection 13.41 provide us with the parameters 
o^expa, (^expcn and Pctra (rccall that actra IS givcu by Eq. (|3.13p ). and therewith with an 
independent way of determining the proportionality constant \/aexpa(3expa + ctctraPctra 
of Eq. 14.81 Repeating this fitting procedure for many impact velocities results in Fig. [T2b . 
where \J o-expaiiexpa + cictra/^ctra IS plotted as a function of log^o Fr. A weak (logarithmic) 
dependence on the Froude number is revealed. It can also be seen that the value Ci « 1.10 
of the proportionality constant in Eq. (|4.8p found from the fit to the closure depth data 
in Fig. [T2]is consistent with the data when one wants to disregard the Fr dependence. 

4.2. Air Entrainment 

After pinch-off, an air bubble is entrapped, as is clearly visible in Figs. [Tb and [2] The 
(rescaled) volume of this bubble VbubUe/hQ is not only found to solely depend on the 
Froude number, but also to exhibit close to power-law scaling behavior: The scaling 
law for the volume of the bubble observed in experiment and simulation is found to be 
Vbubbie/hl cx Fr\ with A = 0.78 (see Fig. [Hi). 

This is surprising, since for the impact of a liquid mass on a free surface the vol- 
ume of air entrained in th e process scales with a different exponent VbubUe oc Fr ' 
( Prosperetti fc Ogua ( 19971 )). In this section we will try to shed light onto the origin of 



this scaling behavior using our findings of Section 14.11 

In Section 14.11 it was found that the axial length of the enclosed bubble at pinch- 
off scales roughly as {zcoii — Zdisk,coii)/ho « 1.10 Fr '", if we ignore the weak Froude 
number dependence of the prefactor ^Jaexpal^expa + oictraPctra (see Fig. [12b V Therefore, 
the scaling of the axial length Zcoii — Zdtsk, coii of the enclosed bubble cannot by itself 
account for the observed scaling of Vtubbie- The radial length scale hrad of the bubble 
must therefore be Froude number dependent and should scale as 



''rad 

OC 



ho 



Vt 



bubble 



•^{)\^ disk, coll ^coll) 



1/2 



CX 



Fr' 



0.78 



Pj,0.50 



1/2 



= Fr^-^* . (4.11) 



Now what would we expect based on our simplified model? The maximum radial 
expansion of the cavity at any depth z is given by h„iax{z), see Eq. (I3.17[) . As the 
depth Zmax at which the radial size of the bubble is maximal is somewhere between 
the closure depth Zcoii and the depth of the disk at closure Zdisk,coii, we have Zmax ~ 
{zcoii + Zdisk,coii)/2 = {i/2)y/aexpaf3expa + ctctraf^ctraFr^^^ ■ If wc insert this depth into 
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hmaxiz), Eq. ()3.17|) . we find 



if^rad CK llmax\Zmax } ^^ i^O \ -L T" , — -TT 

y "^y/ ^expaPexpa i OictraPctra 



/ioVl + 0-26Fr^/^ (4.12) 



In the last (approximate) equation we have used that y/ UexpaPexpa + occtraPctra ~ 1.10 
and oil^pg^Pexpa ~ 0.40 (cf. Fig.[T2b). If a power-law fit hmd/ho vs. Fr is enforced on this 
dependence in the regime 2.5 < Fr < 250 one obtains the observed effective exponent 
0.14, hrad/ho oc Fr"'^ . Alternatively, by taking the square of Eq. (|4.12p and multiplying 
with the vertical extension {zcoii — Zdisk. coii) of the bubble we find the following prediction 
for the bubble volume 

Vbubbie ^hl{l + 0.26Fri/2) Fr^/^ . (4.13) 

Clearly, the model predicts power-law scaling only in the limit of large Froude numbers. 
Moreover, as in this ImatVhubUe oc Fr, the scaling prediction is in agreement with the 
Prosperetti fc Qgua ( 19971 ) result. Again, in the regime 2.5 < Fr < 250 the effective 



exponent is 0.78. 

We test the above prediction by looking at three different quantities that capture the 
radial expansion of the cavity in experiment and numerics. The first is the effective, or 
average, radius hf.ff of the bubble which is computed directly from the experimental and 
numerical cavity profiles (i.e., without any scaling assumption of the axial length scale) 
at the pinch-off time by 

hlg=- — r°" h\z)dz. (4.14) 

\Zdisk. coll ZqqIIJ J Zdisk coll 

The second quantity we look at is the maximal radius of the bubble hjnax.coii ^^ the 
time of pinch-off which is a more direct measure of the expansion of the cavity, hmax. coii 
can be directly observed from the cavity profile at the time of pinch-off as the maximal 
radius for a depth between Zdisk,coii and Zcoii- 

In Figure [T4b we compare these two quantities heff/ho and hmax,coii/ho with a third, 
namely the measured Vbubbie compensated for the expected scaling behavior of its axial 
extension Zdisk,coii - Zcoii oc /loFr^'^, i.e., Vbubbie / hl¥v^ ' "^ . AU of these three quantities 
follow the same trend, which is well described by the prediction Eq. (j4.12p from the 
model (the blue dashed line in Fig. [T4b). and close to the expected Fr*''^^ scaling which 
is indicated by the solid red line. Finally, comparing the measured bubble volume 
Vbubbie I h'^ with the model prediction Eq. (|4.13ll in Fig. [T4k (dashed blue line), we find 
excellent agreement. 



5. Conclusions 

In this article we investigate the purely gravitationally induced collapse of a surface 
cavity created by the controlled impact of a disk on a water surface. We find excellent 
agreement between experiments and boundary integral simulations for the dynamics 
of the interface, as well as for the flow surrounding the cavity. The topology and the 
magnitude of the flow in the simulations agree perfectly with the PIV results. 
In experiments it is found that a secondary air effect, the "surface seal", has a sig- 
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nificant influence on the cavity shape at high Froude number. Since the surface seal 
phenomenon (and its influence) is more pronounced at higher impact velocities, it limits 
our experimental Froude number range. In the boundary integral simulations the air 
was intentionally excluded, thus avoiding this limitation. 

Because the velocity of the impacting disk is a constant control parameter in our 
experiments, a simple theoretical argument based on the collapse of an infinite, hollow 
cylinder describes the key aspects of the transient cavity shape. 

This model accurately reproduces the dynamics of the cavity including its maximal 
expansion and total collapse time. It also captures the scaling for the depth of closure 
and the total depth of the cavity at pinch-off, and predicts their ratio to be close to 2, 
where 2.1 is found in experiments and simulation. 

There is a clo s e sim ilarity of this description to the cavity dynamics proposed by 
Duclaux et all pOOTl l. However, by introducing the asymmetry between the radial 



expansion and collapse, we find a better agreement between the theory and the radial 
dynamics of the cavity. The fact that the flow is qualitatively different during expansion 
of the cavity on the one hand and its contraction on the other is f ound to be responsible 
for the asymmetry. Our approach is also conceptually different, as lDuclaux et al\ ( 20071) 



take a^xpa to be independent of the Froude number, while we allow it to be weakly 
dependent on Froude and, more importantly, our description includes the last stage of 
the collapse, which is solely driven by inertia. 

We find the volume of air entrained by the impact of the disk to behave as 
ybubbie/hQ oc (1 + 0.26Fr^' ^)Fr^' ^ . This dependence is set by the Froude dependence of 
two length scales, namely the axial length scale, distance between the pinch-off point 
and the disk, and the radial expansion of the cavity. Here we have excellent agreement 
between the experimental and numerical findings and the prediction of the model. 

Finally, th e appendices deal with the time evolution of the cavity radius we discussed 



previously in lBergmann et al\ (|2006l ). In this paper, and subsequent papers dealing with 
the universality of the last stag es o f the pinch-off from a theoretical point of view 
(JGordillo &: Perez-SaboridI (|2006l ) and lEggers et all (|2007( )). the pinch-off is assumed to 



be symmetric around the closure depth, whereas in experiment and numerics we observe 
that the minimal radius of the void actually moves downward in time. As this (small) 
axial translation could be relevant for this universality issue, it is studied in detail in 
Appendix |B1 where we find that it can be included within the model presented in this 
paper, as a secondary effect. 
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Appendix A. Revisiting IBergmann et al\ 1 20061 ) 



In t his appendix we will review the results which were presented in iBergmann et ali 



( 20061 ) as far as they are necessary for the understanding of the material in Appendix [Bl 



together with some additional results. 
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A.l. Neck radius at closure depth 

The most prominent length scale describing the cavity dynamics close to pinch-ofF is 
the neck radius. It can be taken either at its minimum (see also Appendix [B]) or at the 
constant depth of closure at each instant of time. In this A ppendix we will deal wi th the 



constant deptn ot closure at eacn mstant oi time, in tins A ppenaix we will deal wi tx 
latter, closely following the discussion in our earlier letter (JBergmann et al\ ( 2006[ )) 



The neck radius hcoii over the time r is found to obey a power law scaling (see Fig.fTSbl. 
where the exponent is observed to vary between 0.55-0.62, depending on the Froude 
number (Fig. 115b). So, for all Froude numbers the scaling exponents are above the value 
of 1/2 that is expected from Eg. (13.81). This result is consistent with the recent careful 
experiments of lThoroddsen et al\ (|2007l ) , for the somewhat different experiment in which 



an air bubble pinches off from an underwater nozzle. 

The deviation can be partly understood by considering the full two-dimensional 
Rayleigh-type equation Eq. ()3.4I) instead of only the fir st term as was d one i n the 
derivation of Eg. (|3.16l) . The procedure is described in iBergmann et al\ ( 20061 ) and 



iGordillo et al\ (|2005l ) and introduces a logarithmic correction to the neck radius, 

Koii{t) ■ {-\og{h,ou)f'^ ^ T^'\ (Al) 

However, even though this result improves the description of the experiments and 
numerics, small deviations are still seen in the dynamics of hcoii{t) for small Froude 
numbers. These deviations, which depend on the Froude number, show the influence of 
the initial conditions on the early stage of the pinch-off. 



As was described in iBergmann et al\ ( 2006 ). these deviations suggest that the neck 



radius is not the only relevant length scale for the cavity shape around the pinch-off 
point. As the cavity shape in axial direction is approximately parabolic, the second 
characteristic length scale can conveniently be chosen as the radius of curvature R in the 
z direction, which is defined as 

1/R{t) = d\ldz\=,_^^^. (A 2) 

We found the time dependence of the radius of curvature to also follow a power-law 
with a Froude-dependent exponent aR. The scaling exponents of the neck radius and 
the radius of curvature differ significantly from one another at small Froude number, 
but tend to converge to 1/2 for higher Fr (see Fig. [T5k and b). 

A. 2. Cavity shape at pinch-off 

If the pinch off would be self-similar, the free-surface profiles near the closure point 
at different times (Fig. 116b ) would superpose when rescaled by any one characteristic 
length, e.g., the neck radius ^.(t). Due to the different time-dependence of the neck 
radius and the radius of curvature, however, such an operation fails to collapse the void 
profiles onto a single shape. To characterize the free-surface shapes, one thus needs to 
consider the neck radius hcoii{t) and the radius of curvature R{t). 



As was argued in IBergmann et al\ (|2006[ ). if the radial dimension r is scaled by hcoii, it 



follows from the locally parabolic shape and Eq. (|A2|1 that the axial dimension z should 
be scaled by ^/hcoiiR- Instead of using the actual length scales to rescale the profiles in 
Fig. 116b . we use the power laws that were obtained from the numerical simulations to 
collapse all experimental profiles onto a single curve. This signals once more the excellent 
agreement between the simulations and the experiment, as the scaling exponents are 
obtained from the simulations and the profiles from experiment. 
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As will be discussed in Appendix [B] the axial position of the minimum radius Zmin 
moves down somewhat as the collapse progresses and it is therefore necessary to translate 
the profiles in the vertical direction to match the position of the mi nimum radius. For 
the in ertially driven pinch off of a bubble from an underwater nozzle 'Thoroddsen et_i 



(120071) also find that two length scal es are necessary to collapse the fr ee su rface profiles 



Rec ent theoretical calculations by iGordillo &: Perez-SaboridI (|2006f ) and lEggers et all 



([20071) for the symmetric inertial pinch off of a single bubble (either initiated by 
surface tension or a straining flow) indicate that the scaling law for hcoii{T) has an 
exponent which slowly varies with time, i.e., strictly speaking it does not scale. Our 
findings so far cannot confirm or disprove this theory, since our experiments and 
boundary integral simulations do not have sufficient temporal range to find the small 
deviations in the exponent. To escape the limitation of the experimental range set by 
the viscosity, surface tension, and air, the experiment should be scaled up to an unre- 
alistic size (with a container size of over 10^ metersljj) to match th e needed prec i sion o f 



at least 10 decades in time presented in the numerical calculations of lEggers et al\ (|2007l ) 



Appendix B. Minimal neck radius 

Recently, a lot of attention has been given to the time evolution of the minimal 
cavity radius hminiT) at the closure d epth and the ( non-) u niversality of its scalin g 
behavior approaching the pinch -off, see Gordillo et ali (120051) . iBergmann et ali (|2006[ ). 



Gordillo &: Perez- Sab or idl (|2006[ ). an dlEggers et all (l2007l) (cf. also A ppen dix E). Espe- 



cially in the theoretical analysis of ICordillo fc Perez-SaboridI (J2006h and lEggers et al 



( 20071 ) it is a key assumption that the dynamics is symmetric around the pinch-off point 



In this Appendix we will show that in our setup this assumption is not satisfied, as the 
position of the minimum neck radius actually moves downwards in time. 

In Figure [TTk the difference 6z{t) between the depth at which the cavity radius is 
minimal Zmin (t) and the depth of closure Zcoii is plotted as a function of the time interval 
T remaining until pinch-off for four different Froude numbers. Each data series is obtained 
from the boundary integral simulations and, as we focus on the behavior close to pinch- 
off, starts when the minimal radius is equal to the disk radius /iq. 

We observe that the dynamics near pinch-off can be reasonably well described by a 
simple proportionality dz oc r. This becomes even more clear if we compensate Sz by 
T in a double logarithmic plot of the same data (Fig. [TTb). which reveals a plateau in 
which Sz/t becomes independent of time, especially for the lower three Froude numbers. 

As the pinch-off is approached, not only the depth of the minimal radius Zmi„(r) 
converges to the depth of closure Zcou, but naturally also the minimal radius itself /imm(T) 
approaches the radius at the depth of closure hcoiiir) as can be seen in Fig. [TTb . The 
relative difference Sh{T) between hmin and hcoii is seen to be smaller for increasing Froude 
number due to the cavity taking a more cylindrical shape at higher Froude number. 

Since the cavity profile is locally parabolic close to the pinch-off point the approach 
of Knin to hcoii is dcscribcd as: 

where R{t) is the radius of curvature in the axial direction at the minimum neck radius. 

t The current experiments, in a 0.5 x 0.5 x 1 m"^ container, cover 2 orders of magnitude in 
/i(t) whereas 10/2 — 5 orders of magnitude are needed. 
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This radius of curvature was found to e xhibit power-law scaling in time with a 
Froude-dependent exponent R{t) oc r"^ (cf. iBergmann et al\ (J2006l ) and Appendix [Aj 
Fig. [T^ which we can combine with the linear time dynamics found for 5z in Fig. [T71 
This leads to hmin approaching hcoii as 



5h oc r 



2-aH 



(B2) 



Indeed, in Fig. [TTH this scaling is confirmed for the lower Froude numbers. For the 
highest Froude number (Fr = 200) the data seem to deviate from this scaling due to 
small deviations from 5z (x t which are observed for this Froude number in Fig. 117b . 



The final question we want to address in this section is whether it is possible to 
understand the relation 5z (xt from the model. We start from Eq. (|3.16p with tcoii = Ai 
given by Eq. (|i^ . i.e., for tcross <t < tcoii 



h{z,t) = .y2hoactraVy/At{z) - t 



(B3) 



To find the depth of minimal radius we now have to compute the derivative of h{z, t) to 
z and equate to zero, or, equivalently 
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dz 
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(B4) 
(B5) 



which leads to the conclusion that the depth of the minimal radius is always equal to the 
closure depth, i.e., independent of time (cf. Eq. (|4.8p l. Clearly this is in disagreement 
with the observations. This was to be expected as the translation of the depth of the 
minimal radius is quite small, typically an order of magnitude smaller than other length 
scales, e.g., the closure depth). We will therefore have to look for a second order effect. 



To take the discussion one step further, we return to Fig. 112b . where we find that 
Oictra slightly decreases with increasing Froude number. From this we can infer (at least 
qualitatively) that adra also very slightly decreases with depth. This stands to reason, 
as Uctra shows the same trend (in Fr) as aexpa , which is the ratio of the initial expansion 
velocity of the cavity and the disk velocity. As at greater depth the hydrostatic pressure 
is larger it is expected that aexpa should decrease with depth. 

Now letting actra (and the other parameters aexpa, Pexpa, and (3ctra) depend on z we 
have 

' dactra 



^[h{z,tf] = 2hQV 



dz 



T + 2hQVactra{z) { ^ 

\ oz 



0. 



(B6) 



For any fixed Froude number Fro we can now Taylor expand dAt/ dz around the closure 
depth Zcoii which leads to 



dAt _ /dAt\ 



d 



(Sz) 



d^At 
dz^ 



(Sz) 



92 At 
dz^ 



(B7) 



as within the model the closure depth is defined by the condition {dAt/dz)z^„ii = (see 
subsection l4.ip . Recall that Sz = z — Zcoii- Combining the above two equations then leads 
to the following relation between Sz and r 



Sz = - 



1 



f daetra 
actra V dz 



d^At 
dz^ 



(B8) 
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where it is good to note that At does not only directly depend on z, but also indirectly, 
through the parameters Qctra, aexpa, Pexpa, and Pctra- 

From the data presented in Fig. 112b . we can perform a local second order fit to adra as 
a function of yPr around \/Ftq. Assuming that the relation between closure depth and 
Froude number also holds near the closure depth, i.e., z/ho = CiVFt with Ci w 1.10 
(see again subsection l4.1[) we can translate this into a quadratic expression in Sz 



Sz (Sz)"^ 

Ctctra ^ actra,Q + ki-. h fc2 ,, ; (B 9) 

ho ftp 

and similarly for the other parameters a^xpai Pexpa: and Pctra- 

As an example we take Frg = 5.1 for which we find actra.o — 0.56, fci = —0.087, 
and fc2 = 0.038. After a straightforward but quite lengthy calculation we evaluated the 
prefactor in Eq. (|B8p to give Sz/Hq « 0.20TV/ho. This is in the same direction and of 
the same order of magnitude as the result from our boundary integral simulation where 
we find Sz/ha w O.blrV/ho (see Fig. [T7k and b). Repeating this procedure for Fr — 10.2 
Fr = 200 yields proportionality constants of 0.16 and 0.074 respectively, which correctly 
predict the downward trend with increasing Froude number that is also observed in the 
numerical simulation, where we find 0.42 and 0.25 respectively. Considering the large 
number of approximations made in this calculation and the subtlety of the effect the 
agreement is remarkable. 
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Figure 13. (a) Double logarithmic plot of the depth at which the void collapses z^ou and the 
depth Zdiak.coii of the disk at collapse time for experiments with four different disk radii (see 
legend) and for the boundary integral simulations (open circles), all as a function of the Froude 
number. Experiments in which a surface seal occurs during the collapse are indicated by red 
symbols, the experiments without a surface seal by blue symbols. Here, a surface seal is said 
to occur if at some point in time due to air suction the splash closes onto itself and the cavity 
is completely sealed off. The experiments without a surface seal are found to agree well with 
the numerically obtained values (open circles) and the theoretical prediction for the scaling 
of Zcoii/hf) = CiFr^'^ (red dashed line) and Zdisk,coH/ho = 2CiFr^'^ (blue solid line), with 
Ci — 1.10 obtained from a fit to the data of Zcoii- The experiments for which a surface seal 
occurs are seen to slightly deviate from this prediction. 

(b) The ratio of the depth of the disk at the time of pinch-off Zdisk, coii and the pinch-off depth 
ZcoU for different disk radii as a function of the Froude number. The experiments without a 
surface seal (blue symbols) agree well with the numerical results (open circles). The ratio for 
the numerical result and experiments without a surface seal lie close to the predicted value of 
2 indicated by the black dashed line. The experiments in which a surface seal occurs are again 
indicated by the red symbols and found to deviate more with increasing Froude number for a 
fixed disk size. The inset shows the definition of the depths ZcoU and Zdisk,coii at the closure 
time. 
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Figure 14. (a) The volume of the bubble VbubUe entrained during the collapse of the cavity from 
experiments (colored symbols) and simulations (black open circles), normalized by the cubed 
disk radius /ig, as a function of the Froude number in a double logarithmic plot. The data suggest 
a power- law scaling VhubUe/ho oc Fr"^ where a linear best fit through the data between Fr — 2.5 
and 250 gives A ~ 0.78 (solid red line, shifted for clarity). The dashed blue line corresponds to 
the model prediction Eq. (|4.13p . 

(b) Double logarithmic plot of three quantities that measure the radial length scale of the 
entrapped air bubble. From bottom to top: the efi'ective (or average) radius h^ff of the bubble 
at pinch-off (Eq. (|4.14p l. the maximum radius of the bubble hmax,coii at pinch-off, and the 
square root of the bubble volume compensated for the expected scaling of its vertical extension 
[Vbubbie/hoFr^'^]^'^ , all compensated with the disk radius ho- The dashed blue line is the model 
prediction Eq. (|4.12p and the solid red line represents a power law with the scaling exponent 
A = 0.14 expected from Eq. (|4.1ip . 
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Figure 15. (a) Time-evolution of the neck radius at closure depth hcoiiitau) and the radius of 
curvature R{t) in the axial direction for Fr = 3.4, 5.1, 10.2, and 200. Both are found to follow 
a power law with respective exponents an and oir. (b) Power-law exponents an and a_R of the 
neck radius and radius of curvature as a function of the Froude number. Here we have added 
the exponents ay which are corrected for the logarithmic factor in Eq. (|A l[l . In the large Fr 
limit, all exponents approach the Rayleigh- value of 1/2 [cf. Eq. (|3.16|) ]. 
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Figure 16. The profiles of the void obtained by experiments for ho = 30 mm and V — 1.0 m/s 
(Fr = 3.4). (a) Void profile obtained at different instances in time, (b) Void profile in which 
the radial and axial coordinates were rescaled with the powerlaws of hcoii (t) and \/hcoiiit)R{t) 
respectively. Here, the numerically determined power-laws were used for the neck radius hcoii (t) 
and the radius of curvature R{t) (see main text). 
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Figure 17. (a) The difference (5z(r) between the depth of the minimal radius Zminij) and 
the depth of closure ZcoU (both normalized by the disk radius fto) as a function of the time 
interval t (normalized by fto/V^) remaining until pinch-off for different Fronde numbers. Each 
data series is obtained from the boundary integral simulations and starts when the minimal 
radius equaled the disk radius /iq. (b) Doubly logarithmic plot of the same data as (a), only 
now &z is compensated by r revealing an interval in which &z is proportional to r (showing 
as a horizontal line in the plot), (c) The difference 5h{T) = hcouij) — hmin(T) of the radius 
at closure depth hcoiiir) and the minimal radius hmin(T) of the cavity, normalized by hcou(T), 
as a function of r (normalized by ho/V) for different Fronde numbers. Clearly, Sh/hcoii is 
smaller for increasing Froude number which results from the fact that the cavity shape becomes 
more cylindrical for high Fr. (d) Doubly logarithmic plot of the same data as in (c) , only now 
hcoii — hmin is Compensated by r^ to reveal the resultant scaling r~°''^ (see text). The dashed red 
lines are the scaling exponents a_R(Fr) found by determining the time evolution of the radius 
of curvature in the axial direction through the procedure outlined in iBergmann et al] (|2006l ) 
(afl(Fr=3.4) = 0.29; afl(Fr=5.1) = 0.32; aH(Fr=10.2) = 0.33; and aH(Fr=200) = 0.40)). The 
jumps in the data of Figures (a) and (c) are a result of the regridding routine in our boundary 
integral simulation and have no physical meaning. 
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